MPP-2011-11 



The meV mass frontier of axion physics 

Georg G. RafFelt,^ Javier Redondo,^ and Nicolas Viaux Maira^ 

^ Max- Planck- Institut fiir Physik (Wemer-Heisenberg-Institut), Fohringer Ring 6, 80805 Munchen, Germany 
^ Departamento de Astronomia y Astrofisica, Pontificia Universidad Catolica de Chile, 
Av. Vicuna Mackenna 4860, 782-0436 Macul, Santiago, Chile. 
(Dated: 19 August 2011) 

We explore consequences of the idea that the cooling speed of white dwarfs can be interpreted in 
terms of axion emission. In this case the Yukawa coupling to electrons has to be gae ^ 1 x 10~^^, 
corresponding to an axion mass of a few meV. Axions then provide only a small fraction of the cosmic 
cold dark matter, whereas core-collapse supernovae release a large fraction of their energy in the form 
of axions. We estimate the diffuse supernova axion background (DSAB) in the universe, consisting 
of 30 MeV-range axions with a radiation density comparable to the extra-galactic background light. 
The DSAB would be challenging to detect. However, axions with white-dwarf inspired parameters 
can be accessible in a next generation axion helioscope. 
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I. INTRODUCTION 



II. COOLING OF COMPACT STARS 



The Peccei-Quinn (PQ) mechanism remains perhaps 
the most compelling explanation for the absence of CP- 
violating effects from the QCD vacuum structure [H-Q- 
An unavoidable consequence is the existence of the ax- 
ion, the Nambu-Goldstone boson of a new U(1)pq sym- 
metry. Axions acquire a mass rua ^ 'mT^f-,^/ fa by their 
mixing with neutral mesons, where niTr = 135 MeV and 
/tt = 92 MeV are the pion mass and decay constant, 
and fa is a large energy scale related to the sponta- 
neous breaking of U(1)pq. Axions generically interact 
with hadrons and photons. They may also interact with 
charged leptons, the DFSZ model [3] being a generic case. 
All interactions are suppressed by , so for large /q, ax- 
ions are both very light and very weakly interacting. Re- 
actor and beam-dump experiments require < 30 keV 
0, while precision cosmology excludes the rUa range 
1 eV-300 keV Sub-eV mass axions would still be 
copiously produced in stars. The cooling of white dwarfs 
(WDs), neutron stars and supernova (SN) 1987A pushes 
the limits to nia < 10 meV [l], i.e. fa > 10^ GeV. Wc 
here explore the impact of axions near this limit, i.e. the 
mcV frontier of axion physics. 

This range is complementary to the other extreme of 
the allowed axion window. During the QCD epoch of 
the early universe, the axion field gets coherently ex- 
cited, generating a cold dark matter (CDM) fraction of 
Pa/PCDM e?(10 peY/maY-^ 0, where 9; = ai//a 
is the initial "misalignment angle" relative to the CP- 
conserving value. For Qi ^ 1, axions with rUa ^ 10 ficY 
{fa ^ 10^^ GeV) provide all of CDM and can be detected 
in the ADMX experiment Q. If the reheating temper- 
ature after inflation was large enough to restore the PQ 
symmetry, our visible universe emerges from many do- 
mains and an average (Qf) ^ tt^/B has to be used. In 
this case, axions also emerge from the decay of topolog- 
ical defects and the CDM density could correspond to 
nia as large as a few 100 fieV @. Either way, meV-mass 
axions provide only a subdominant CDM component. 



The most restrictive astrophysical limits on those ax- 
ion models that couple to charged leptons arises from 
WDs. An early study used the WD cooling speed, as 
manifested in their luminosity function, to derive a limit 
on the axion-electron coupling of gae ^ 4 x 10~^^ (loj . 
In the early 1990s it became possible to test the cool- 
ing speed of pulsating WDs, the class of ZZ Ceti stars, 
by their measured period decrease P/P. In particular, 
the star G117-B15A was cooling too fast, an effect that 
could be attributed to axion losses if ^ 2 x 10"^'^ [11|. 
Over the past twenty years, observations and theory have 
improved and the G117-B15A cooling speed still favors 
a new energy-loss channel [l^l- What is more, the WD 
luminosity function also fits better with axion cooling if 
= 0.6-1.7 X 10-13 [g. 

While complete confidence in this intriguing interpre- 
tation is certainly premature (perhaps even in the need 
for a novel WD cooling itself), the required axion pa- 
rameters are very specific, motivating us to explore other 
consequences based on the WD benchmark. 

Axion cooling of SNe has been widely discussed in the 
context of SN 1987A d, The 10 s duration of 

the neutrino burst supports the current picture of core 
collapse and cooling by quasi-thermal neutrino emission 
from the neutrino sphere. New particles that are more 
weakly interacting than neutrinos, such as the axions dis- 
cussed here, can be produced in the inner SN core, leave 
unimpeded, and in this way drain energy more efficiently 
than neutrinos, which can escape only by diffusion. The 
SN 1987A neutrino burst duration precludes a dominant 
role for axions. Quantitatively, this argument depends on 
the model-dependent axion-nucleon couplings, the uncer- 
tain emission rate from a dense nuclear medium, and on 
sparse data. As we shall see, the limit does not preclude 
the WD interpretation, but a SN would lose a significant 
fraction of its energy in the form of axions. 

The speed of neutron-star cooling as measured by the 
surface temperature of several pulsars [l7j is another pos- 
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sible laboratory to search for axion coohng and a hmit 
comparable to the SN 1987A bound was found How- 
ever, neutron-star cooling depends even more dramati- 
cally on nuclear physics uncertainties and on the details 
of axion-nucleon coupling than the SN 1987A bound so 
that it is hard to make such arguments precise. How- 
ever, if the WD interpretation applies, axion emission is 
another effect to be taken into account in the complicated 
theory of neutron-star cooling. 



III. DIFFUSE SN AXION BACKGROUND 

Returning to the energy loss of proto neutron stars af- 
ter core collapse, axions saturating the SN 1987A limit 
are emitted as copiously as neutrinos. Then one not only 
expects a strong axion burst from each SN, but also a 
large cosmic diffuse background flux from all past SNe, 
the diffuse SN axion background (DSAB) in analogy to 
the diffuse SN neutrino background (DSNB) [H. All 
past SNe in the universe provide a local flux of order 
fO cm^^ s^^ [l^ that will become detectable in a Gd- 
enriched version of Super-Kamiokande [1^ or a future 
large scintillator detector [2l| with a rate of a few events 
per year. The estimated core-collapse rate is scaled to the 
amount of extra-galactic background light (EBL), repre- 
senting the integrated star-formation history [22|. The 



intensity of the EBL is 50-100 nW m~^ ster~^, corre- 
sponding to an energy density of 13-26 meV cm~^, i.e. 
about 10% of the energy density provided by the cosmic 
microwave background. 

The present-day average core-collapse rate is i?cc = 
1.25 X 10-4 Mpc"^ yr-i and increases with redshift 
roughly proportional to 10^ until z = 1 and then flat- 
tens or slightly decreases [l^- Assuming that every SN 
releases 3 x 10^'^ erg in the form of neutrinos of all flavors 
and integrating over i?cc, properly redshifting the energy, 
leads to a present-day DSNB of 26 meV cm"'^, almost 
identical with the EBL. In other words, stellar popula- 
tions release on average as much gravitational binding 
energy in the form of neutrinos as they release nuclear 
binding energy in the form of photons. 

For meV-mass axions, therefore, the energy density of 
the DSAB can be comparable to the DSNB and the EBL, 
and indeed would be the most important axion popula- 
tion in the universe. The axion losses of ordinary stars 
would contribute a much smaller energy density, just as 
the neutrinos emitted by all ordinary stars contribute an 
energy density of only about 7% of the EBL [H, [H. 

The DSAB will be calculated in analogy to the DSNB 
where the Ve spectrum as a function of present-day De 
energy E is given by the redshift integral 



dN 
dE 



dz 



{{l + z)^m + z)]}{R,,iz) ^ }, (1) 



where i?cc is the core-collapse rate. The function ip{E') 
provides the number of De per rest-frame energy inter- 
val dE' released by an average SN. Notice that A'tot = 




CO [MeV] 

FIG. 1. Diffuse backgrounds of SN neutrinos and axions, 
assuming that either one carries away the full SN energy 
(3 X 10^^ erg). The width of the bands reflects only the uncer- 
tainty in the core-collapse rate Rcc- For a thermal spectrum 
with T = 4 MeV is assumed, carrying away 1/6 of the total 
energy, whereas for axions we use the bremsstrahlung-inspired 
spectrum of Eq. (|8j with Tcoro = 30 MeV. 



/ dE (1 -I- z)if[E{l + z)], the total number of neutrinos 
released by a SN, is invariant against redshift. Further, 
\dt/dz\-^ = Ha{l + z)[nA + fluil + zf]'^^^ with cos- 
mological parameters taken as Hq = 70 km s~^ Mpc~^, 
Ha = 0.7, and JIm = 0.3. Notice that \dt/dz\ and i?cc ac- 
tually form one combined factor proportional to the ratio 
of the average luminosity per galaxy in SN neutrinos rel- 
ative to stellar photons. Assuming Etot = 3 x 10^^ erg, 
1/6 of this in De, and Tp^ = 4 MeV after flavor oscil- 
lations [11], we show the^DSNB in Fig. [H The width 
of the band reflects only the uncertainty of Rcc, not the 
uncertainty of SN neutrino emission. 



IV. AXION PROPERTIES 

Before estimating the DSAB we summarize the rele- 
vant phenomenological axion properties. Their mass is 

™ - ^ "^-^^^6mevH!^, (2) 
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fa 



fa 



where z = mu/md = 0.35-0.60 [3|, but we always use the 
canonical value z = 0.56. The interaction with fermion 
/ has the axial- vector derivative structure 



(3) 



where Cf is a numerical coefficient and gaf = Cfirif/fa 
the corresponding Yukawa coupling. For protons Cp = 
[Cu - 1/(1 + z)]Au + [Cd - z/{l + z)]Ad and neutrons 
Cn = [Cu - 1/(1 + z)]Ad + [Cd - z/{l + z)]Au, where 
Au = 0.84 ±0.02 and Ad ^ -0.43 ±0.02 Q. The axion- 
photon interaction is 



r 



(4) 
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with = a/(27r/,) [E/N - 2/3(4 + z)/(l + z)] « 
al{2iTfa) {E/N - 2) and E/N is the ratio of the elec- 
tromagnetic and color anomalies. 

We assume DFSZ axions 3 for which E/N = 8/3, 
C„ = i sin^ /? and Cd = Ce = ^ cos^ ^. The WD value 
gae = 1 X 10^^'^ for the axion-electron coupling implies 



fa = 1.7 cos^/? 10^ GeV. 



(5) 



We choose cos^ /3 = 1/2 because much smaller values 
would favor SN emission over WD emission and lead to 
overly optimistic DSAB estimations. This choice implies 
fa = 0.85 X 10^ GcV, ma = 7 meV, Ce = 1/6, C„ = 



and Co 



-1/3. Then ga^ 



1.0 X 10"^^ GeV"' is 
two orders below the CAST sensitivity [2^. The axion- 
nucleon couplings are gan = and gap = 3.7 x 10~^". 

V. ESTIMATING THE DSAB 

For such small couplings, axions escape freely from a 
SN core once produced. The dominant production pro- 
cess is nucleon bremsstrahlung NN — )■ NNa, but a re- 
liable calculation has proven elusive @. Axions couple 
to the nucleon spin and are produced in spin fluctua- 
tions caused by the tensor force in NN collisions. Early 
calculations used nondegenerate free nucleons and a one- 
pion exchange (OPE) potential in Born approximation. 
However, if the nucleon spin-fluctuation rate T„ were as 
large as found in these calculations, destructive interfer- 
ence from multiple scattering would reduce the emission 
rate [lB|. On the other hand, based on measured NN 
scattering data, the spin-flip cross section was found to 
be much smaller than implied by the OPE approxima- 
tion [li]. 

For a phenomenological description of the axion inter- 
action with a nuclear medium we use the dynamical spin- 
density structure function S'cr (a;) p^[l6j. The absorption 
rate for axions of energy Ea is g^pf {p/8m^) uj Sa{Ea) 
where p is the matter density. If we ignore spin cor- 
relations between different nucleons, the normalization 
is J/^^ Sa{i^) duj/ {2Tr) = 1. Emission and absorption 
are related by detailed balancing, implying S^{—uj) ~ 
5cr(a;) e~"/^. The spectral axion emission per unit vol- 
ume is therefore 



dEa " 16 7r2 ml^"^ ' 



(6) 



The energy- loss rate per unit mass Q/p is given by 
ta = (SaTv/^^) C^^/"^?/) ^ in. terms of the dimensionless 
integral F = Sa{~Lo) iijj /Tf dw /2tx . 

Low-energy bremsstrahlung is essentially a classical 
phenomenon. Classical spins kicked by a random force 
with a rate F^ imply Sa = Fcr/(w^ -I- r^/4) and inspire a 
one-parameter representation fulfilling all requirements 



5*0.(0;) 



+ r2/4 'Z^ + l 



(7) 



We consider only interactions with protons (typical abun- 
dance 30% per baryon) and adopt F = 1 as a rough esti- 
mate so thatea - 1.6x10^7 crgg^^ s-i(T/30 MeV)^. 
The SN 1987A neutrino signal duration requires ta % 
1 X 10^^ erg g~^ s~^, providing gap < 0.8 x 10~^ 

Our assumptions correspond to Fo- ~ 2T . Assuming 
an isothermal SN core, the axion spectrum is 



El 



dn 



dEa El + r2 eEJT + I 



(8) 



This distribution is of course not very well determined 
and mostly serves the purpose of illustration. With 
Tcorc ^ 30 MeV we find {Ea) - 80 MeV. 

For our WD inspired axion parameters, a SN emits 
roughly 1/8 of its energy as axions. Considering all un- 
certainties, this fraction could be smaller or as large as 
1/2, at which point it would seriously affect the SN 1987A 
signal. Assuming all energy is emitted in axions and with 
the spectral shape of Eq. we find the DSAB shown in 
Fig. [T] The average axion energy is about 35 MeV. 



VI. DETECTING THE DSAB? 

Detecting this flux is extremely challenging. Axions 
with the parameters considered here interact much more 
weakly than neutrinos of comparable energy. The DSNB 
will be detectable in Super-Kamiokande and in next gen- 
eration large-scale detectors, but the DSAB produces a 
much smaller signal. 

One may think that conversion in large-scale astro- 
physical magnetic fields may provide a detectable signal. 
It would have to stick above the diffuse gamma-ray back- 
ground in the 30 MeV region that was measured by the 
EGRET satellite to be 1-2 x 10"^ cm^^ g-i sr"i MeV"^ 
[131. Dividing the DSAB in Fig. [T] by 47r to obtain a flux 
per sterad, we see that the conversion probability would 
have to be of order 10~*. 

In a transverse B field and after travelling a distance 
L, the axion-photon oscillation probability is 



Pa^^ = {ga^B/qf sm\qL/2), 



(9) 



where q = {rn^ — m'^)/2E. For nia = 7 meV we can ne- 
glect the photon plasma mass in interstellar space. For 
E = 30 MeV the oscillation length inE/ml is 1500 km. 
For these parameters, the maximum conversion rate is 
Pa^f = 6 X 10"^^ (B/Gauss)^, apparently too small 
for any realistic astrophysical B-field configuration. For 
axion-like particles (ALPs), in contrast, where nia and 
ga-y are independent parameters, large conversions and 
astrophysical signatures are conceivable [2^, [l^ . 

The situation does not improve for a — >■ 7 conversions 
near compact objects such as pulsars or active galactic 
nuclei where B fields can be much larger and the pho- 
ton plasma mass can be such that q = js^]. The 
Hillas diagram of possible sources of high-energy cos- 
mic rays shows that BxL^lOGxpc can be attained. 
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The maximum conversion probability (taking g — >■ 0) is 
{ga-yBL/2)^ ~ 0{1). However, the intrinsic 7-ray emis- 
sion tends to be far too large to disentangle the two com- 
ponents, even if spectral features could help (sTj . 



VII. NEXT GALACTIC SN 

The next galactic SN will provide a high-statistics 
signal of 10 MeV range neutrinos [13 ■ What about 
the comparable energy release in 100 MeV range ax- 
ions? The DeP — >■ ne~^ cross section is cFp^p ^ 9.4 x 
10"''^ cm^ (i?p^/MeV)^, so for detection energies of 20- 
30 MeV it is around lO"**" cm^. For axions, a reaction 
like a+p N + tt has a cross section on the A resonance 
{Ea - 340 MeV) of order 100 mb {U/faf ^ IQ'^^ cm^. 
Most of the axion flux is not on resonance and the rates 
for such reactions always seem too small for realistic de- 
tection. The largest conceivable SN signal is in a future 
megaton detector and if the red supergiant Betelgeuse 
at a distance of 200 pc collapses. This scenario provides 
about 3 X 10^ Ve events and conceivably a few events 
above the tail of the neutrino spectrum that could be at- 
tributed to ap — > Nit. Of course, such an scenario would 
require a much more careful discussion. 

VIII. NEXT GENERATION AXION 
HELIOSCOPE 

A more realistic detection possibility of WD-inspired 
axions is with a large helioscope beyond CAST. The con- 
version probability is (5a^i?i/2)^ 10^^^ for L = 20 m, 
B = 10 T and = 10" GeV"^ The solar axion 
flux from processes involving electrons is 0.47 x 10~^Lo 
with an average energy of 2.1 keV and a flux at Earth of 
2.0 X 10^ cm~^ s~^ 133|, yielding several events per year 
and m^. The feasibility of such an instrument with an 
aperture up to 4 m^ has been recently assessed [33 |. 

One amusing application for such an instrument is to 



detect axions from a possible Betelgeuse SN explosion. 
Assuming all SN energy is released in axions of average 
energy 80 MeV, Betelgeuse provides an axion fluence at 
Earth of 5 x lO^'' cm~^. In this case one needs an aperture 
exceeding 20 m^ to get a few events. Pointing the instru- 
ment at Betelgeuse in time for the explosion is possible 
by the early warning few days) provided by the de- 
tectable thermal neutrinos from the silicon burning phase 
preceding core collapse [ssj . 

IX. CONCLUSIONS 

The intriguing hint from WD cooling for the exis- 
tence of DFSZ-typc axions with fa ^ 10^ GeV and 
ma 7 meV implies that core-collapse SNe emit a large 
fraction of their energy as axions. The universe would 
be filled with 30 MeV-range axion radiation with a den- 
sity comparable to the diffuse SN neutrino background 
and the extra-galactic background light. The axion pop- 
ulation produced in the early universe would comprise 
only a small fraction of cold dark matter, but of course 
the cold dark matter in the universe may well consist 
of different components. Searching for a sub-dominant 
meV-mass axion dark matter component is a new chal- 
lenge that has not yet been seriously addressed in the 
literature. It is intriguing that axions with such parame- 
ters are accessible in a next generation axion helioscope, 
a possibility that should be vigorously pursued. The in- 
terpretation of WD cooling in terms of axion emission is, 
of course, speculative, but it suggests a fascinating new 
mcV-mass frontier of axion physics. 
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